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Abstract 

Consider a simple two-way relaying channel where two single-antenna sources exchange information via a 
multiple-antenna relay. To such a scenario, all the existing works which can achieve full diversity order are based on 
the antenna/relay selection, where the difficulty to design the beamforming lies in the fact that a single beamformer 
needs to serve two destinations. In this paper, we propose a new full-diversity beamforming scheme which ensures 
that the relay signals are coherently combined at both destinations. Both analytical and numerical results are 
provided to demonstrate that this proposed scheme can outperform the existing one based on the antenna selection. 



Index Terms 



Two-way relay systems, beamforming, network coding, amplify-and-forward, diversity, outage probability. 
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I. Introduction 

Relaying technique has been considered as an efficient method to extend the coverage and improve the 
system throughput. Particularly two-way relay systems have drawn a lot of attention, since it can double 
the spectral efficiency of information exchanging by applying network coding [[T)-[|3]. By broadcasting the 
mixture of both source messages at the relay, the communication between two sources can be accomplished 
in two time slots, where each source recovers its desired message by subtracting self-interference from 
its observation. Compared with the traditional 4 time-slot schemes, network coding can increase the 
throughput significantly. 

To further improve the transmission performance, multiple -input and multiple-output (MIMO) tech- 
niques have been introduced to two-way relay channels By properly designing the beamform- 
ing/precoding, multiple data streams can be transmitted simultaneously in MIMO systems, which provides 
another reliable approach to improve the performance of relay transmissions. Specifically the performance 
analysis of network coding in two-way amplified-and-forward (AF) relay systems was derived in [@]. A 
joint source and relay precoding design for two-way AF relay systems was studied in JH, where the 
multiplexing gain can be maximized by channel parallelization. In (61, an interference alignment scheme 
is proposed for multiple two-way relay systems, and the inter-stream interference is eliminated by applying 
precoding, while the intra-stream interference is coped with by using network coding. 

Although there have been such extensive studies about two-way relaying channels, some challenging 
issues are still left as open problems. For example, consider a simple two-way relaying channel where 
two single- antenna sources exchange information via a multiple-antenna relay. To such a scenario, all 
the existing works which can achieve full-diversity gain are based on the relay/antenna selection, and 
to the best knowledge of the authors, how to design a full diversity relay beamformer is still an open 
problem. The motivation to design the beamforming is analogue to the performance gain of maximum ratio 
combining (MRC) over the selection detection in single-input multiple-output (SIMO) systems. However, 
in the context of two-way relaying, the difficulty to design the beamforming lies in the fact that a single 
relay beamformer needs to serve both destinations. In this paper, we focus on a simple two-way relaying 
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scenario and the main contribution of this paper can be summarized as follows: 

Firstly, a full-diversity relay beamformer is designed for the addressed scenario, to which the conven- 
tional beamforming/precoding techniques, such as the generalized singular value decomposition (GSVD)- 
based precoding in 0, cannot provide any type of spatial transmission gain. The key idea of our proposed 
scheme is to utilize the symmetry of the observation phase at two destinations. As a result, the relay signals 
can be coherently combined at both destinations. Just like the performance gain of MRC over the selective 
detection, the proposed beamforming scheme can achieve the same diversity gain as the existed antenna 
selection scheme in 0, Q, and also offer better reception reliability in the low signal-to-noise ratio 
(SNR) region. 

Secondly, to evaluate the performance of the proposed transmission scheme, the outage probability 
is analyzed in this paper. Particularly an upper bound of SNR is first developed, which facilitates the 
development of a closed-form expression of outage probability. The asymptotic analysis of such an 
bound is also provided to demonstrate the full-diversity gain achieved by the proposed scheme. Finally, 
the simulation of outage probability shows that our proposed scheme outperforms the antenna selection 
scheme, especially when the SNR is low. In addition, the numerical results also demonstrate the accuracy 
of the derived analytical results. 

The rest of this paper is organized as follows. Section II describes the system model, and introduces the 
proposed full-diversity transmission scheme. In Section III, the performance analysis of proposed scheme 
is studied. The simulation results are shown in Section IV, and followed by the conclusions in Section V. 

Notation: Vectors are denoted as boldface small letters, i.e., a, a m denotes the m-th element of a and 
a T is the transpose of a. |c| is the Frobenius norm of c, and c can be either a vector or a number. E{x} 
is the expectation of the variable x, and Pr{a; < p] denotes the probability that the value of the random 
variable x is less than p. T(d) is the Gamma function, 7(0,2;) denotes the lower incomplete Gamma 
function, K u (z) is the modified Bessel function with imaginary argument, and all the referred special 
functions follows the form given in ffTOTl . 
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II. System Model and Protocol Description 

Consider a two-way relay system, where two sources Si and S 2 exchange messages via a relay R. As 
illustrated in Fig. \T\ each source node is equipped with one antenna, while the relay is equipped with 
N antennas. All nodes are assumed to employ the time division duplexing for simplicity. Due to the 
symmetry of time division duplex systems, the incoming channel and the corresponding outgoing channel 
are symmetric. All the channels are modeled as quasi-static Rayleigh fading channels, and the relay has 
the access to full source-relay channel statement information (CSI). 

The transmission can be accomplished in two phases by applying network coding. During the first phase, 
both sources transmit their own messages to the relay simultaneously, and the network coded observation 
at the relay can be expressed as 

y R = hx t + gx 2 + n R , (1) 

where X\ is the message transmitted by S\, and its transmit power is limited as E{xf x{\ = P\, h denotes 
the Af x 1 channel vector for the link from S x to R, x 2 and g are defined similarly for S 2 , and the 
corresponding transmit power is P 2 , and n R is the additive Gaussian noise at the relay. Assuming that the 
AF strategy is applied, the relay broadcasts its network coding message after beamforming in the second 
phase, and the observations at Si and S 2 can be expressed as follows respectively, 

Vl = /3h T Qhx! + (3h T Qgx 2 + (3h T Qn R + m , (2) 

V2 = /3g r Qgx 2 + /3g T Qhxx + f3g T Qn R + n 2 , (3) 

where Q is the beamforming matrix at the relay, (5 = ■^ |Qhpp 1 +|c^|2p 2+ |Q|2 .2 is the power normalization 
factor at the relay, P R is the relay transmit power and rij is the additive Gaussian noise at Si, i = 1,2. 
By using the global CSI, Si can subtract its self-interference from the observed network coding message, 
and obtains the desired message from the other source. 

A. Full-Diversity Beamforming Design for Studied Two-Way AF Relay Systems 

Since the relay is with multiple antennas, antenna selection is a straightforward method to achieve 
full diversity gain, and improves the transmission performance. On the other hand, when all the relay 
antennas are utilized, the conventional beamforming/precoding schemes, such as the generalized singular 
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value decomposition (GSVD) based scheme in Q, cannot achieve the best diversity performance, which 
is due to the fact that the signal space at the sources is flattened into one dimension. To provide a full- 
diversity gain scheme with better performance, the beamforming design is proposed in this subsection. 
First the definition of diversity is introduced as follows, 

Definition 1: Spacial diversity can be achieved by redundantly receiving the same message over inde- 
pendent fading channels. Particularly the diversity gain is a if the transmission scheme satisfies that 

pout 

lim — — = a (4) 

where P ont describes the probability that the transmission fails, as defined in [8J, and p denotes the average 
SNR. 

To have a closer look at the diversity order, we further derive the definition of outage probability based 
on the results given in [|S), 

pout = pr{i? < Rth} = Pr{SNR < 7th}; (5) 

where R is the transmission data rate, i? t h is set as the threshold for R, SNR is the receive SNR and 7 tn 
is its threshold. The last equation shows that both the outage probability and the diversity gain are mainly 
determined by the receive SNR. Recalling the observations given in © and ©, the SNRs for S\ and S2 
can be derived as follows respectively, 

(3 2 \h , vol . ,„ , vol 

SNR 1 = T^!9 , -, N 9 = nl^l 19 , !/-> 19 , 1^19 / Pi ( 6 ) 



/3 2 |h T Qg| 2 P 




|h T Qg| 


2 


(/3 2 |h T Q| 2 + l)a 2 


2|Qh| 2 - 


f |Qg| 2 


+ IQIVp 


/3 2 |g T Qh| 2 P 




|g T Qh| 


2 


(/3 2 |g T Q| 2 + l)a 2 


2|Qg| 2 - 


f |Qh| 2 


+ IQI7p 



SNR2 (/3 2 |g T Q| 2 + l)<r 2 2|Qg| 2 + |Qh| 2 + |Q| 2 /p P ' (7) 
where the transmit power is set as Pr = P\ — Pi — P for simplicity, and the average SNR p = P/a 2 . 

Limited by the single antenna setting at the source nodes, the key step to achieve full diversity gain 
is to design the relay beamforming matrix Q. As introduced in [9|, the SNRs of full-diversity schemes 
share the common structure as follows. Over ai^xl channel with the fading vector f, 

K 

SNR ful - div = |f|V = X)lMV (8) 

k=l 

Such a result shows that the diversity gains can be achieved by coherently accumulating the signal power 
from different pathes at the destination, which coincides with the clarification in Definition 1. Then we 
focus on the signal parts in the SNRs of S\ and S%, 

N N N N 

h T Qg = J2 Yl <l™nh m gn = Yl (lmn\h m \\gn\e Mm+en \ (9) 

m=l n=l m=l n=l 
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N N N N 

g T Qh =J2Y1 IrnnQmK = ?mnbm| I h n\ S 9 ™*^ , (10) 

m=l n=l m=l n=l 

where <\> v denotes the argument of h p and 6 q is defined similarly for g g . As shown in © and (flOl) . 
the uncertainty of arguments blocks the reception enhancement at each source node. To ensure the 
beamforming design can benefit both sources, the beamforming matrix can be defined as 

r e ~^^ m+6n \ m = n 

Q = {<lrnn)NxN = \ ■ (11) 

0, m 7^ n 

Substituting ([IT]) into © and ©, the SNRs at Si and 5*2 can be further derived as 

CATD ( Em=l \hm\ \9m\) oattd ( Sm=l \hm\ \9m\) 

bNRi = r- — ——. — p. bNR 2 = , lo — -771- P- (12) 

2|h| 2 + |g| 2 + N/p F) 2|g| 2 + |h| 2 + N/p y 

By adjusting the argument of channel fading for each path, the receive signal power can be strengthen 

at both sources, and thus the full-diversity gain can be achieved, which is demonstrated in the next section. 

III. Performance Analysis for the Proposed Full-Diversity Transmission Scheme 

In this section, the performance of the proposed full-diversity transmission scheme is evaluated. Firstly, 
a tractable upper bound of outage probability is derived, and then its closed-form expression can be 
obtained. Next the asymptotical high SNR approximation of the derived upper bound is analyzed, which 
demonstrates the diversity gain of the proposed transmission scheme. 

A. Upper Bound of Outage Probability for the Proposed Scheme 

Recalling the definition given in ©, the outage probability can be derived by applying the probability 
density function (PDF) of SNR. Since the channel vectors for both links are coupled tightly in the signal 
parts of SNRi and SNR 2 , it is hard to derive their PDFs. To find a feasible method to evaluate the outage 
probability of the proposed transmission scheme, we first derive a tractable upper bound for the SNRs, 
and the following lemma is presented. 

Lemma 1: Denoting that w = — -, — (^ m=1 l fem l l gm ^ can ^ e b OU nded as follows. When the 

number of relay antennas iV is large enough, 

A^oo 

w > 1. (13) 

Proof: We first further derive w as 

_ |(^m=i l^-mPl^ml 2 + 2 ^2 m ^ n \h m \\g m \\h n \\g n \) 

w ~ TL v N \h |2T7-Lv^ i 12) ' 

\N 2^m=l \ n ™\ )\N 2^m=l \9m\ ) 



IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY, VOL. XX, NO. X, X XXXX 



As introduced previously, h m and g n are independent Gaussian distributed variables. Based on the law of 
large numbers, the denominator of w in (fl4)) can be approximated as follows, 



(jt E i^i 2 ) (I E \^\ 2 ) N =°° NM 2 ) + *J Kki 2 ) + *J = (i + w + $o. d5) 

^ m=l ' ^ n=l ' 

where 5i and 5 2 approaches zero when iV goes to infinity. Similarly, both parts of the numerator in (fl"4l) 
can be given as follows when N approaches infinity, 

1 N 

^ \ h m\ 2 \9m\ 2 ^Ed^H^I 2 ) +5 3 = l + 5 3 , (16) 

m=l 

^Yl l^m||^||^n||^| Af = 00 2(iV-l)E(|/ im ||^ m ||/ in ||^|) +5 4 = i(iV-l)7r 2 + 5 4 , (17) 

where the last equation follows the fact that \h p \ and \g q \ are independent Rayleigh distribution variables, 
and 5 3 and <5 4 are infinitesimals. Substituting (fl"5T ). (fT6l and (fTTT) into (fl4l) . we can obtain that 

1 + ()i + 2 + 0lt>2 ^ 

The lemma has been proved. ■ 
Based on Lemma 1, the lower bound of the receive SNRs can be given as follows. When the number 
of relay antennas iV is large, the channel vectors can be decoupled at the signal parts for both sources, 



SNR 1 > ar/mi.., , ll , at/ n P» SNR 2 > Ar/ol , ' , , ,' , iT77 ^ - d9) 





|h| 2 |g 


2 


iV(2|h| 


2 + 


g| 


2 + iV/p) 





|hp 


|g 


2 


iV(2|g| 


2 + 


h| 


2 + iV/p) 



Then an upper bound of outage probability can be derived in closed-form, which is presented in the next 
theorem. 

Theorem 2: The upper bound of outage probability for the proposed scheme in Section II can be 
expressed as 

7(iV,^) ^/iv_ A(2A r Tth )i. 



r(N) 3 J v{N) 



2Na 

1 V I I I I ( 

N-j, 



(N-j)b e 



i=0 k=0 1=0 

where a = and b = Nj t h(l + 27th)- 



3=0 

N-l i N-l /.\ /Ar 1X A rt? ; N+k-l 4 < f _ k /n v _( 2JV+1 >Tth 
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Proof: Without loss of the generality, the upper bound of outage probability for S\ is derived. Denoting 
that x\ = |h| 2 and x 2 = |g| 2 , it can be given as 

pout-up _ p i 3:1X2 ^ 0/ 

P - Fl \N(2x 1 + x 2 + N/p) P<lth 

= Pr {x 2 < 2Na] + Pr <j < x x < lth%2 * f a , x 2 > 2Na\ = K x + K 2 . (21) 

As introduced previously, x\ and x 2 are two independent Chi-squared distribution variables, whose PDFs 
can be shown as 

-xi ~N—1 P —X2 

f = ^1 - f = ^ ! (22) 

J XI p^ ,JX2 p^Q • K^) 

Then substituting (1221) into (121) . /Ci and fc 2 can be obtained, 

^e- 7(iV,2iVa) 

£2 = / /m ( / * ) dx 2 



2 

2A^a WO 



r(jv) 

1 f°° ( N 7t h x 2 + Na 



2 I A^ ZZoAr'S. K-V^cfafr. (24) 



7 th x 2 +A f a 



To obtain the closed-form expression, (|24l can be expanded as follows by denoting that z 



px 2 -2N lt ^ 



' -e p z -Tth dz — 



nmJa ( P z- lth ) N+1 

b ^ [°° z\2N lth z + Na) N - 1 (2N^ X+Na \ 
1 ^ iV J i=0 A T(i){pz -7th) 

where /C 2 is simplified by applying 7(^,2) = T(N)[l - e- z (J2iU f%)] in ED- Then we focus on 

deriving the closed-form expression of C 1 and £ 2 respectively. Particularly denoting that t = pz — 7^, 

£x can be derived as 

k -2A-. p gjv 7lhf + 6)^-1 » ^ /AT -IV e- M « /~ e"* 

= fwl e * = i j (2JV7,h)6 Wl i^n* 



AT-l 



v ,iv-i\ (2iy 7th y 
2s \ n I r(N\ 

j=0 



j J r(iv) 



iV-j. 



-2Na 

—-, (26) 
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] dt 



] dt 



where C\ is simplified by using Binomial theorem. Following the notations in (1261 ), the closed-form 
expression of C 2 can be obtained in a similar way, 

AT 1 (2JV+l)7 tv , 

^ - r(iv)r(z)p^ J 

N-l % N-l , . s , N , . , jv+fc-i i+i _ fc , x - (2JV+1)7 th 

Then the proof has been finished. ■ 
Although the derivation of p out - u P is based on the assumption that the number of relay antennas is 
large, the simulation results in the next section show that the derived upper bound is a general expression, 
i.e., p out - u P is also applicable when N = 2. 

B. Analysis of Diversity Gain for the Proposed Scheme 

To derive the achievable diversity gain of the proposed scheme, the asymptotical analysis of p out ~ u P 
is given. Specifically it can be obtained by applying high SNR approximation on the original expression 
of p out ~ u P given in (fUT) . and the derivation is given in detail as follows, 

1) High SNR Approximation of ' K\ in (1211) : For high SNR region, 2N J th approaches 0, and thus the 
lower incomplete gamma function in /Ci can achieve the following asymptotic form [fTOl . 



Then the high SNR approximation of /Ci can be derived as 



(2AT 7th )^ 1 
NT(N) p l 



= \J; b J, — ■ (29) 



2) High SNR Approximation of /C 2 in (I2TT) .- To derive an accurate approximation, we begin with the 
expression of /C 2 given in (f24l) . When p is in the high SNR region, the argument of the lower incomplete 
Gamma function in (1241) approaches 

7th^2 + (Nq/th/p) = lthX 2 + (iVWp) _^ 7th (3Q) 
px 2 - 2iV 7th p{x 2 - (2jV 7th /p)) p ' 
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Then the approximation of /C 2 can be derived as follows by substituting (|30b into (|24l ). 



Similar to the approximation of K\, the asymptotic form of /C 2 can be given as follows, 



l-7|iV,^ 

P 



(3D 



(2Nj 



th) 



< 1 . J 1 



+ o Tv ■ (32) 



Based on the high SNR approximation of /Ci and /C 2 , the asymptotic form of p out ~ u P can be given in 
the following corollary, which demonstrates the diversity gain of the proposed scheme. 

Corollary 3: For the high SNR region, the derived upper bound of outage probability for the proposed 
scheme, which is provided in Theorem 2, can be approximated as 

NT(N) p N N[T(N)] 2 p N \p N J ' ^ J 

Substituting (|33l) into ©, it is easy to show that the diversity gain of the proposed scheme is N, which 
is the full-diversity for the studied scenario. 

IV. Numerical Results 

In this section, the numerical results are provided to evaluate the transmission performance and demon- 
strate the accuracy of theoretical analysis in the previous section. In Fig. [21 the cumulative distribution 
functions (CDF) of w defined in Lemma 1 are plotted when the number of relay antennas is set as 
N = 3,4, 5, 6. As shown in the figure, Pi{w > 1} raises as N increases, and approaches 1 when N ^ 3. 
Such results verify the asymptotic inequality given in Lemma 1. 

Fig. [3] gives the numerical results for the derived upper bound of outage probability and the correspond- 
ing asymptotic analysis, where N is set as 2, 3 and 4. The simulation results show that the derived upper 
bound p out ~ u P is quite tight in the high SNR region, especially when N is small. Moreover, the outage 
probability of the proposed scheme decreases faster when the relay is equipped with more antennas, which 
demonstrates the diversity gain. The curves for the upper bound and the asymptotic analysis have the same 
slope as the Monte-Carlo curve, and thus the analysis results in Section III can be testified. 

To further evaluate the outage performance of the proposed scheme, the comparison with the antenna 
selection scheme and the direct AF scheme for the studied scenario is provided in Fig. @] and the numbers 
of relay antennas are given as N = 2, 3, 4, respectively. Since the direct AF scheme transmits the messages 
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without additional operations, it does not obtain any diversity gain, while both the proposed scheme and 
the antenna selection scheme can achieve full-diversity gain. As shown in the figure, the proposed scheme 
can achieve better performance than the antenna selection scheme, especially in the low SNR region. 
Compared with the antenna selection scheme, our proposed scheme can reduce the required SNR by 
1.3 dB with N = 3 and when the outage probability achieves 1CT 2 , and such performance gain can be 
enlarged as iV increases. 

V. Conclusion 

To achieve full-diversity gain, a joint beamforming and network coding scheme was proposed for two- 
way relay systems. By combining the messages from different paths at the sources, the transmission 
reliability can be improved. Then a closed-form upper bound of outage probability was derived for the 
proposed scheme, and the high SNR asymptotic analysis was also given to demonstrate the achieved 
diversity gain. The simulation results were provided to verify the derived analysis results, which also 
shows that our proposed scheme outperforms antenna selection scheme in the studied scenario. 
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Fig. 1. System model. 
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Fig. 2. The CDF of w in Lemma 1. 
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Fig. 3. The derived upper bound of outage probability, TV = 2, 3, 4. 
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Fig. 4. The outage probability performance, N — 2,3, 4. 



